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Введение. Актуальность исследования контактных задач 
для цилиндрических тел обусловлена необходимостью 
проведения инженерных расчетов на контактную проч- 
ность валов, стержней и трубопроводов. В настоящей ра- 
боте изучается новая контактная задача статической тео- 
рии упругости о взаимодействии жесткого бандажа с бес- 
конечным двухслойным цилиндром, состоящим из внут- 
реннего сплошного и внешнего полого цилиндров, между 
которыми выполняются условия гладкого контакта. 
Наружный цилиндрический бандаж посажен с натягом и 
имеет конечную длину. При помощи интегрального пре- 
образования Фурье задача сводится к интегральному 
уравнению относительно неизвестного контактного дав- 
ления. 

Материалы и методы. Рассматриваются разные комбина- 
ции линейно-упругих материалов составного цилиндра. 
Исследуется асимптотика функции-символа ядра инте- 
грального уравнения в нуле и бесконечности, играющая 
важную роль для использования аналитических методов 
решения. Для решения интегрального уравнения вводится 
основной безразмерный геометрический параметр и при- 
меняется сингулярный асимптотический метод. 
Результаты исследования. В соответствии со свойствами 
функции-символа предложена специальная легко факто- 
ризуемая аппроксимация этой функции, пригодная в ши- 
роком диапазоне изменения параметров задачи. При по- 
мощи метода Монте-Карло рассчитаны параметры этой 
аппроксимации. Получены асимптотические формулы как 
для контактных давлений, так и для их интегральной ха- 
рактеристики. Расчеты сделаны для разных материалов и 
относительных толщин цилиндрического слоя, в том чис- 
ле для тонкостенных слоев. 

Обсуждение и заключения. Полученные асимптотические 
решения эффективны для относительно широких банда- 
жей, когда размер области контакта превышает диаметр 
составного цилиндра. Важно, что используемый метод 
остается применимым и для случаев, когда внешний ци- 
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линдрический слой можно рассматривать как цилиндри- 
ческую оболочку. Асимптотические решения можно ре- 
комендовать инженерам для анализа контактной прочно- 
сти упругих деталей цилиндрической формы с упругим 
покрытием из другого материала. 
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составной цилиндр, аппроксимация, асимптотика. 
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Введение. Динамическая контактная задача для преднапряженного упругого цилиндра, наполненного 


жидкостью, изучалась в работе [1]. Статические контактные задачи для однородных упругих тел цилиндриче- 
ской формы рассматривались в работах [2—6] при помощи регулярного и сингулярного асимптотических мето- 
дов. Было установлено [4], что для цилиндрических тел символы ядер интегральных уравнений контактных 
задач характеризуются более сложным асимптотическим поведением в нуле и бесконечности, чем, например, в 
контактных задачах для упругой полосы. Это потребовало применения усложненных аппроксимаций этих сим- 
волов легко факторизуемыми функциями при использовании сингулярного асимптотического метода. Особен- 
но усложняется аппроксимация для полых тонкостенных цилиндров [6]. Предложенная аппроксимация [6] эф- 
фективна даже для случаев, когда тонкостенный упругий цилиндр можно рассматривать как цилиндрическую 
оболочку [7]. Исследовалась контактная задача о взаимодействии упругого кольца с упругим цилиндром [8]. 
Износ упругого цилиндра анализировался в работе [9]. Цель настоящего исследования — получить решение 
контактной задачи для составного двухслойного упругого цилиндра на основе сингулярного асимптотического 


метода и эффективной аппроксимации символа ядра интегрального уравнения. 

Материалы и методы. В цилиндрических координатах и, 2 (при осевой симметрии) рассмотрим беско- 
нечный упругий составной цилиндр внешнего радиуса А, который состоит из внутреннего сплошного цилиндра 
радиуса А!<А с упругими параметрами у1, С, (коэффициент Пуассона и модуль сдвига) и внешнего цилиндри- 
ческого слоя с упругими параметрами у, С. Между слоем и внутренним цилиндром выполняются условия сколь- 
зящей заделки. Рассмотрим контактную задачу о взаимодействии описанного составного цилиндра с жестким бан- 
дажом по области |2<а. При заданном натяге бандажа 6 требуется определить контактные давления 

с, =-9(2) (25а). 

Используя интегральное преобразование Фурье для решения собственно смешанной (контактной) краевой 

задачи для уравнений Ламе упругого равновесия и вводя безразмерные обозначения (штрихи далее опускаем) 








Ат, 5' -5, Н -=, ик, = к я (1) 
получим следующее интегральное уравнение относительно 4(б): 
ДЕ =лб (|(<1), (= [дшдсовбидаи, (2) 
где символ ядра имеет вид ” | 
п --Ч чак, (8) 
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—А,5.455 [А (Аз Ава В Аз Аз ) =. Ар (АзАв = А Ав У, 


4, =—(А, Аи - А.А, А (АА, -А,Ак)- ААА] + 
+А, А [А (Аз А, -А,А:)-А (А.А, - АА], 

А =и\ -2(-у)Л, А, ==—иК. -2(-у)К, А, =и|, А. =—иК, 
А, = (3—2у)ш — (и? +4 -у))1, А, =—(3-2\иК, — (м? +4А-у))К,, 
Аз =и(-и), Ад =—и(К! +иК,), 

А; =2(1-у)Г", А, =2-У)К!, А, =—4А-ум Г, 
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А =иКГ —=21-у)", А’ =-иКК, -2А-У)К'", Аз =ий, Ад 
А; =и К =-2-\,), Ак =иш, 
А, = 2 и -Е+4А-УЕГ, 
А, =—-—2иК: — (и +4а-УЕУКУ, 
Аз =ик Ч и? , Ад =-иКИ К! и? КЕ, 
А, = 6—2)" —Рк+4-У Е, А = ви -шЬ), 
=Г(и), К, =К, (и), Г =Г (К), К” =К, (иК), п=0, 1. 


Здесь //(и), К‚(и) — модифицированные функции Бесселя. Безразмерный параметр ^, характеризует от- 
носительную ширину области контакта. 
Функция Г(и) в нуле и бесконечности ведет себя следующим образом: 


(У-1@+=-\, +5%,) +”, 
2(у-1[+10+=-у, +=%,)+А"=] | 
& = (У+0(\, -9-=(-—П(\, +1. (4) 





Низ (м) = 240) = 


гр 
Ци) = —+—+0(и”) (и—> +0), Р=1-2у. 
ии 


При А=0 значение Г(0) совпадает с известным для однородного сплошного цилиндра [4]. 

Для решения уравнения (2) применим сингулярный асимптотический метод [3,4], эффективный при 
достаточно малых значениях /.. С целью применения метода Винера-Хопфа [10] использовалась легко факто- 
ризуемая аппроксимация функции Г(и) (3) выражением 


ми? + В? Бр и? + А*С? В ур) 
: 5_ехр| == ‚ @ = ехр| — |. 
и +С ) [2 +10“ 10 











-. +? 10) — 


При расчетах брали два случая: 1) железо внутри, цинк снаружи (&=2,126, у=0,27, у,=0,28); 2) 
алюминий внутри, цинк снаружи (&=0,779, у=0,27, у,=0,34). В таблице 1 даны значения параметров 
аппроксимации (5), ее относительная погрешность на действительной оси 09 (%), рассчитанные при 
использовании метода Монте-Карло при разных относительных толщинах внешнего слоя К. 

Таблица 1 
Параметры аппроксимации 





А В С 0 А В С 0 


К Железо внутри цинка Алюминий внутри цинка 


0,10 1,230 1,549 4,638 2,5 0,451 6,113 3,394 
0,30 1,291 1,563 4,042 2,5 0,888 5,977 3,029 
0,50 1,296 1,851 5,006 3,0 2,399 1,467 2,524 
0,70 1,258 2,057 6,376 3,0 0,720 7,945 2,512 
0,90 1,024 4,360 92,381 5,0 1,481 9,890 9,369 
0,99 3,540 24,178 5,347 7,0 1,178 5,535 | 235,946 10 
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Поскольку, как следует из (4), 
1-у%, 
э 

21 -—у)(1+\,) 
то в случае малых # аппроксимация, рассчитываемая по формуле (5), для тонких внешних слоев должна услож- 
няться путем увеличения числа входящих в нее параметров. 

Результаты исследования. В результате применения метода Винера-Хопфа главный член 
асимптотического решения интегрального уравнения (2) при малых /^. можно построить в форме 


_6| [1+6 Е 1 
4() т т |4 >. т ((х|51), (7) 


_ #5) + 5) 


И’ ($) = аи, р ег (/В5 +[5-1 О(АС,5) , 


О(Е,5) = Е ехр(-Ез)ег /(В — Е)5) + т ег (Ву). 


Здесь егКх) — интеграл вероятностей. 
Для интегральной характеристики решения 


Р= [ час (8) 


вит (0) = (6) 





Ця) = ори Эк, 0?Эат, 
ИЯ 





на основании формул (7) получим выражение 


Р 2 2 2 2 р: ь 
С 2(2)+= (2 о О эк. (0? эат, (9) 
2 т =) а ехр( вы Е И Сиб, 


Т(Е,5) = < тм =) а ехр(-— во - 
С \ехр(-Р 5) — 
( = ЕЕ еРСВ-Е)5). 


Как показывают расчеты, погрешность асимптотик (7), (9) при ^.<1 не превышает (5+0)%, где 0 — по- 
грешность аппроксимации (5). 
В таблице 2 приведены значения интегральной характеристики Рб’', рассчитанные по формулам (9) 

















при разных значениях Ки ^. 
Таблица 2 


Значения Рб" 


























= 2 1 0,5 0,25 2 1 0,5 0,25 
К Железо внутри цинка Алюминий внутри цинка 
0,10 3,25 5,77 190,9 21,1 3,22 5,71 10,8 20,9 
0,30 3,37 6,02 11,4 22,2 3,20 5,65 19,7 20,7 
0,50 3,67 6,65 12,7 24,8 3,09 5,51 10,5 20,4 
0,70 4,22 7,78 15,0 29,4 2,96 5,33 19,2 19,9 
0,90 5,60 10,3 19,6 38,2 2,93 5,28 9,99 19,4 
0,99 6,91 12,3 23,2 44,9 2,99 5,31 9,94 19,2 



































Заключение. Как ВИДНО ИЗ таблицы яя с уменьшением А, интегральная характеристика контактных 
давлений возрастает, что связано с увеличением площади области контакта. Для случая более прочного 
материала внутри цинка (железа) контактные давления больше, чем для алюминия внутри цинка. С утончением 
слоя цинка вокруг железа (при возрастании 19) контактные давления существенно возрастают. При утончении 
слоя цинка вокруг алюминия этого не наблюдается, поскольку модуль продольной упругости (а также модуль 
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сдвига) у алюминия немного меньше, чем у цинка. Найденные асимптотики можно рекомендовать инженерам 
для анализа контактных прочностных характеристик деталей цилиндрической формы с покрытием. 
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